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SECTION  I 
INTRODUCTION 


A.  GOALS 

Surface  waves  often  arrive  at  stations  located  at  teleseismic 
distances  from  the  epicenter  with  amplitudes  at  or  below  the  noise  level. 

In  order  to  detect  such  signals  and  determine  their  magnitudes,  special 
processing  techniques  must  be  employed.  Three  such  techniques  are 
evaluated  in  this  report  using  data  from  the  Very  Long  Period  Experiment. 
They  are:  chirp  matched  filters  (CMF),  reference  waveform  matched  filters 
(RWMF),  and  the  three- component  adaptive  processor  (TCA). 

The  specific  goals  of  this  study  are: 

• To  estimate  potential  signal-to-noise  ratio  gains  of  each  of 
these  techniques. 

• To  evaluate  the  effectiveness  of  these  three  techniques  in 
increasing  the  surface-wave  detection  capability  of  the  VLPE 
stations. 

• To  apply  the  signal-to-noise  ratio  improvement  estimates  to 
the  calculation  of  surface-wave  magnitudes  for  ev  mts  which 
were  not  detected  on  the  bandpass-filtered  tracr  , 

• To  compare  the  relative  effectiveness  of  the  three  techniques. 


B.  RESULTS  FROM  PREVIOUS  STUDIES 

A preliminary  evahiation  of  these  techniques  as  applied  to 
VLPE  data  was  presented  in  Special  Report  No.  14  (Lambert  et  al.  , 1973), 
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where  the  application  of  each  technique  to  events  having  epicenters  in  a 
small  region  of  Sinkiang  Province,  China  was  discussed.  Conclusions 
were  limited  by  the  small  amount  of  observational  data  available.  The 
preliminary  evaluation  did  indicate  that  appropriate  chirp  matched  filters 
performed  essentially  the  same  as  reference  waveform  matched  filters  when 
matched  with  Rayleigh  waves.  The  use  of  matched  filters  decreased  the  number 
of  non-detected  events  by  36  percent.  Both  chirp  and  reference  waveform 
matched  filters  yielded  signal-to-noise  ratio  improvements  of  3.  5 dB  for 
earthquakes  and  3.  7-3.8  dB  for  presumed  explosions.  The  three-component 

W-. 

adaptive  processor  yielded  detection  results  comparable  to  those  for  the 
chirp  matched  filter. 

Analysis  of  data  recorded  at  the  Alaskan  Long  Period  Array 
(ALPA)  (Strauss,  1973)  indicated  that  chirp  matched  filters  were  slightly 
more  effective  than  reference  waveform  matched  filters,  that  matched 
filters  reduced  the  number  of  non-detected  events  by  20  percent,  and  that 
the  greatest  change  in  the  detection  versus  bodywave  magnitude  plots  caused 
by  inclusion  of  these  detections  occurred  at  the  50  percent  detection  level. 

Analysis  of  data  recorded  at  the  Norwegian  Seismic  Array 
(NORSAR)  (Laun  et  al.  , 1973)  indicated  that  reference  waveform  matched 
filters  were  slightly  more  effective  than  chirp  matched  filters,  that  matched 
filters  reduced  the  number  of  non-detected  events  by  about  10  percent,  and 
that  inclusion  of  these  detections  in  the  detection  versus  bodywave  magnitude 
plots  decreased  the  detection  levels  between  the  30  percent  and  80  percent 
detection  levels.  (The  three-component  adaptive  processor  was  not  applied 
to  ALPA  or  NORSAR  data.  ) 
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c. 


SIGNAL-PLUS-NOISE-TO-NOISE  RATIO  IMPROVEMENTS 


One  of  the  methods  of  comparing  the  performance  of  the  three 
data  enhancement  techniques  under  consideration  was  in  terms  of  signal-to- 
noise  ratio  in  provement  over  the  equivalent  bandpass  filter  (0.023-0.059  Hz) 
signal-to-noise  ratio.  The  signal  - to -noise  ratio  improvement  of  a matched 
filtered  trace  over  the  corresponding  bandpass  filtered  trace,  expressed  in 
decibels,  is: 

IMPROVEMENT  (dB)  = 20  log  (S  /N  )/(S  /N  ) 

10  L'  M M DP  BP'J 


Or,  in  a more  convenient  computational  form: 


IMPROVEMENT  (dD)  = 20  log 


M 


N 


10 


BP 


+ 20  log 


BP 


10  N 


M 


where  S and  N are  the  peak  signal  and  the  RMS  noise  amplitude,  respectively, 
the  M subscript  denotes  matched  filter,  and  the  BP  subscript  denotes  band- 
pass filter. 


Since  the  manner  in  which  the  VLPE  data  was  edited  often 

precluded  the  existence  of  a noise  sample  of  suitable  length  on  the  edited 

signal  trace,  the  values  of  N and  N in  the  above  equation  were  determined 

BP  M 

from  a noise  sample  of  the  same  day. 


Since  the  signals  are  not  noise  fi  ee,  the  signal  amplitudes 
are  actually  signal -plus-noise  amplitudes.  Therefore,  we  will  refer  to  the 
signal-plus-noise-to-noise  ratio  (SNNR)  from  this  point  on. 

In  Section  II  of  this  report,  the  data  base  is  presented  in 
detail,  with  a complete  listing  of  event  parameters,  description  of  available 
stations,  and  discussion  of  the  data  quality.  Section  III  presents  the 
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analysis  of  tlie  cliir))  matched  filter  data,  Section  IV  the  analysis  of  the 
reference  waveform  matclied  filter  data,  and  Section  V of  the  analysis  of 
tlie  three-component  adaptive  processor  data.  Section  VI  presents  the 
comparison  of  the  tlirc'c  techniques  and  conclusions  beisecl  on  these 
c ompa  ri  sons , 
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SECTION  II 


THE  DATA  BASE 

A.  THE  AREAS  OF  INTEREST 

Three  regions  were  initially  selected  for  this  evaluation  of 
the  CMF,  RWMF,  and  TCA  processor  data  enhancement  techniques.  I’hese 
regions  are  the  central  Asia  region  (CENA),  the  Greece- Turkey  region 
(GTUR),  and  the  eastern  Kazakh  test  region  (EKAZ).  These  are  sliown  in 
Figure  II-l.  The  first  two  regions  were  selected  to  evaluate  the  applica- 
tion of  the  data  enhancement  technques  to  earthquakes,  w'ith  the  intent  of 
investigating  any  differences  in  behavior  of  each  technique  due  to  source 
and  path  differences.  Tlie  third  region  (EKAZ)  was  selected  to  evaluate 
the  application  of  these  techniques  to  presumed  explosions.  Unfortunately, 
this  region  yielded  only  ten  presumed  explosions,  which  is  too  small  a data 
base  to  allow  meaningful  analysis.  Therefore,  eastern  Kazakh  presumed 
explosions  are  not  included  in  this  study.  In  future  work,  when  sufficient 
data  is  available,  events  from  this  region  will  be  investigated  in  terms  of 
matched  filters  and  the  three-component  adaptive  processor. 

B.  AVAILABLE  DATA 

A total  of  53  earthquakes  from  central  Asia  and  28  earthquakes 
from  Greece- Turkey  were  processed.  The  parameters  for  these  events  are 
listed  in  Tables  II- 1 and  II-2,  respectively.  The  event  numbers  (ENVO)  are 
taken  from  Lambert  et  al.  , 1973.  All  events  processed  lie  in  the  inmiod 
1 January  1973  to  15  March  1973.  This  time  frame  was  chosen  for  the 
following  reasons; 

II-l 
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• It  contained  sufficient  events  from  the  selected  regions 
to  allow  a meaningful  analysis. 

• A large  number  of  stations  were  operational  during 
this  period. 

• The  amount  of  data  to  be  processed  was  not  so  great  as 
to  preclude  completion  in  the  time  allotted  to  the  task. 

The  stations  used  in  this  analysis  were  CHG  (2),  EIL  (!3), 

KON  (6),  KIP  (8).  ALQ  (9).  ZLP  (10),  and  MAT  (11).  The  se  are  shown  on 
Figure  II-l.  Operational  stations  are  indicated  by  an  X in  Tables  II-l  and 
II-2.  The  remaining  stations  were  not  included  because  they  either  were 
not  operational  during  the  time  period  chosen  or  failed  to  record  any  events 
suitable  for  use  as  reference  waveform  matched  filters.  (It  was  desired 
that  the  entire  data  base  be  processed  by  each  of  the  three  techniques.  This 
could  not  be  done  if  a reference  waveform  matched  filter  was  not  aval’ able 
for  each  station. ) 

In  order  to  evaluate  the  dB  SNNR  improvements  of  each  data 
enhancement  technique,  each  event  as  recorded  at  a given  station  was 
considered  to  be  an  individual  event  and  is  hereafter  referred  to  as  a 
station-event.  The  reason  for  this  designation  is  that  the  travel  paths 
and  epicentral  distances  to  each  station  are  independent  of  each  other. 

Thus,  a chirp  or  reference  waveform  matched  filter  suitable  for  an  event 
recorded  at  one  station  will  not  be  of  value  for  that  same  event  recorded 
at  a second  station.  This  manner  o*"  considering  the  data  yields  a data 
base  of  212  station-events  froni  <"ENA  and  101  station-events  from  GTUR. 

In  order  to  make  comparisons  with  data  reported  previously, 
the  detection  levels  were  computed  in  terms  of  a network,  using  events 
rather  than  station-events.  The  configuration  of  available  stations  conforms 
most  closely  with  those  used  for  Network  3 (Lambert  et  al.  , 1973).  The  two 
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seismic  regions  were  combined  for  this  phase  of  the  study  to  form  a subset 
of  the  data  base  used  to  evaluate  Network  3.  To  decrease  the  possibl’it^  of 
false  alarms,  an  event  was  considered  to  be  detected  only  if  it  was  detected 
at  two  or  more  stations.  By  requiring  two  or  more  stations  to  be  opera- 
tional, a data  base  of  79  events  was  available  for  detection  level  computations. 

The  bodywave  magnitude  distribution  and  maximum-;.ikelihood 
detection  capability  curve  (Ringdal,  1974)  are  given  by  Figure  II-2  for  this 
combined  region.  This  figure  shows  that  for  this  region  and  the  network, 

ju 

the  50  percent  bandpass  detection  level  is  4.72  and  the  90  percent  bandpass 

detection  level  is  5.  25  m units.  We  note  that  Network  3 had  c 50  percent 

L 

bandpass  del  'ction  level  of  4.  50  and  a 90  percent  bandpass  detection  level 

of  4.  99  m units  for  tlie  case  where  detections  on  at  least  two  stations  were 
b 

required.  The  difference  in  detection  level  between  the  two  networks  is 
attributed  to  the  fact  that  the  data  for  Network  3 was  derived  from  a much 
larger  region  than  was  the  data  for  the  network  under  consideration. 

C.  DATA  QUALITY  AND  DETECTION  CRITERIA 

All  station-events  initially  proposed  for  analysis  by  the  three 
data  enhancement  techniques  were  first  screened  for  mixing  and  instrument 
malfunctions.  All  mixed  station-events  were  rejected.  If  only  one  component 
of  a station-event  contained  a malfunction,  it  was  accepted  as  part  of  the 
data  base,  since  it  could  still  be  used  in  the  matched  filter  detection  studies. 

If  more  than  one  component  contained  a malfunction,  the  station-event  was 

Throughout  this  report,  two  decimal  places  will  be  used  when  discussing 
detection  levels  and  M meansurements  in  order  to  maintain  the  convention 
of  earlier  reports.  It  is  recognized  that  the  degree  of  accuracy  in  such 
measurements  actually  permits  only  one  decimal  place. 
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DETECTION  PROBRBILITY 


40 

»— 

z 

UJ 

>■ 

UJ 


u.. 

D 

OC 

Ui 

OQ 

SI 

D 


8-^ 


6-^ 


2-J 


eJ7 


Vm 


/ 

/ 

// 

1///// 


O 


// 


p 

//I 


/ 

1 


r 

yu 


n 


mu 


n T7  1 rr 

3.0  4.0  3.0  ' 

MflGNITUDE  CMB) 
MR  niSTRlBUTION  OF  PROCESSED  EVENTS 


rri  I 


D 

0 


onfcriED 
NOT  DETECTED 


MRGNITUDE  CMB) 


MRXIMLJM-LIKELIHaDD  DETECTRBILITY  CURVE 

FIGURE  11-2 

DETECTION  STATIST! ::S  FOR  THE  COMBINED  REGION 


rejected.  The  results  of  this  screening  process  are  givc-n  in  Tables  II-l 
and  II -Z,  where  an  X undi-r  a station  number  indicates  that  particular 
station-event  was  included  in  the  analysis. 

Lambert  el  al.  , (197  3)  reported  that  erratic  static  gains  were 
encountered  from  lime  to  lime  on  the  horizontal  components  at  \irlually  all 
stations.  Therefore,  SNNR  irnpro ven’iei.ts  and  surface-wave  magnitude 
values  were  calculated  only  for  the  vertical  comjionc'nt  data.  The  horizontal 
data  was,  however,  included  in  the  detection  statistic  s. 

The-  crileri.a  for  determining  whether  a delc'ction  was  achiesu'd 
with  the  use  of  any  of  the  three  enhancement  technicpies  of  a station-event 
a re: 

• A peak  appeared  on  the  matched  filter  or  TCA  processor 
respon  se  trace  3 dB  above  any  other  peak  in  the  first  half  of 
the  signal  gate  defined  by  the  estimated  arrival  time  and  the 
Rayleigh  duration  time. 

• A station-event  was  listed  as  a detection  if  it  was  detected 
on  either  the  vertical  (Rayleigh  wave)  or  transverse  (Love 
v'ave)  components. 

• An  event  was  considered  to  be  detected  if  it  was  detected  at 
two  or  more  stations. 
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CHIRP  MATCHED  FILTER  EVALUATION 


A.  DISCUSSION 

Linear  chirp  matched  filters  were  applied  to  all  the  station  - 
events  listed  in  Tables  II-l  and  II-2.  The  chirp  filters  were  specified  and 
applied  in  the  fr  equency  domain,  using  a chirp  bandpass  of  0,023  to  0,059 
Hz,  Chirp  filter  application  was  performed  by  Fourier  transforming  the  data 
trace  to  the  frequency  domain  and  performing  cross-correlation.  After 
application  of  the  filters,  the  data  were  inverse-transformed  to  obtain  time- 
domain  chirp  filter  outputs. 

The  chirp  filter  response  function  is: 


G (K) 


I: 


i27r(C/N)  (K-K  ) 
o 


if  K < K<  K„ 

Li  H 

ifO<K<K^  orK„<K<N/2 
Li  H 


G (-K)  = G (K)=!' 


where: 


K 

K and  K 

JLJ  H 

K 

o 

N 

C 

G(K)=’^ 


the  discrete  Fourier  transform  frequency  index, 
the  lowest  and  highest  frequencies  in  the  passband, 
the  frequency  index  at  which  zero  phase  shift  occurs, 
the  number  of  transform  points, 

a parameter  which  controls  the  length  of  the  corres- 
ponding time-domain  waveform, 

complex  conjugate  of  G(K) 
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This  yields  a dispersive  time-domain  chirp  vvavc-form  with  a 

linear  group  delay  and  essentially  flat  amplitude  at  all  periods  in  tlie  band 

corresponding  to  K K < K (Harley,  1971). 

Lj  In 

To  minimize  computer  and  Calcomp  plotter  usage,  it  was 
necessary  to  obtain  some  idea  of  optimum  chirp  length  in  terms  of  event 
epicenter- station  separation.  Therefore,  tin*  following  procedure  was 
carried  out  for  the  two  seismic  regions. 

First,  a subset  of  the  data  base  for  each  region  was  formed, 

c:onsisting  of  all  station -events  detected  on  the  bandpass -filtered  trace 

Chirp  filters,  rangl.’.g  in  length  froi  50  seconds  to  1250  seconds  for  ('vents 
o 

closer  tlian  60  , and  from  300  seconds  to  1500  seconds  for  invents  farther 
o 

tlian  60  , v.'ith  length  increment  of  50  seconds,  were  applied  to  each  station- 
event  in  the  subset.  The  lin^e  trace  having  thehighest  peak  value  of  each 
sequence  of  chirp  filler  rt'sponaes,  as  picked  from  the  computer  print-out, 
was  then  plotted  to  ensure  that  it  was  in  fact  a detection.  Plots  of  chirp 
length  versus  epicentral  distance  were  ♦h-'n  made,  using  the  length  associated 
with  the  highest  peak  value.  A least-mean-square  fit  was  then  made  to  the 
data.  This  produced  an  estimate  of  the  (optimum  chirp  H'ngth  in  terms  of 
e|)icentral  distance. 

Next,  the  chirp  length  increments  were  set  as  listed  in 
lable  lll-l.  The  reason  for  changing  the  length  increment  with  increasing 
epii  ( ntral  distance  is  that  if  a constant  increment  were  used,  it  could  form 
too  course  a "grid"  for  close  events,  with  the  result  that  the  o))timum  chirp 
length  would  be  missed.  For  distant  events,  this  same  increment  could  form 
to(j  fine  a "grid,  " resulting  in  wasted  computer  and  plotter  time.  The  chosen 
< hirp  length  increments  were  selected  in  an  effort  to  minimize  this  problem. 
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TADLE  m-1 

CHIRP  FILTER  BOUNDS  AND  INCREMENTS 


Epi  eentral 
Distance- 

Ine  remint 

Lower  Bound'*' 

Uj^i^er  Bound'*' 

< 25"^ 

25  see 

-150  see 

+ 1 50  set 

O C) 

25  -40 

40  sc'c: 

-1  50  sec 

+ 1 50  set: 

o , o 

40  -60 

50  sec 

- 200  s<  c 

4 200  set 

/ o o 

60  -80 

7 5 sec 

- jOO  sec 

4 300  sec 

o 

> CO 

100  sec 

- 300  sec 

4 400  set 

The*  bounds  ar<‘  given  in  senuids  below  (miniis)  or  secends 
above  (])lus)  the  1 east-mean- square  fit  to  the  data. 


Finally,  upper  and  lower  limits  in  the  range  of  chirp 


lengths  were  selected.  These  are  listed  in  Table  III-l. 

The  bounds  were  set  so  that  at  least  85  percent  of  the 
optimum  chirp  lengths  of  the  LR-V  data  (as  determined  from  the  detected- 
on-bandpass  data  subsets)  lay  between  them. 

Using  the  above  information,  chirp  matched  filters  were 
next  applied  to  those  station -events  which  were  not  detected  on  the  bandpass 
filtered  trace.  This  w;ts  performed  in  a two-pass  operation.  For  the  first 
pass,  a suite  of  chirps  was  generated  using  the  appropriate  increment  and 
set  of  bounds.  The  highest  peak  amplitude  of  the  chirp  responses  in  the 
appropriate  time  gate  was  picked  from  the  computer  print-out.  For  the 
second  pass,  the  chirp  length  corresponding  to  this  peak  was  used  as  the 
center  length  of  a set  of  three  chirps,  again  separated  in  length  by  the 
appropriate  time  increment.  The  chirp  responses  so  obtained  were  then 
plotted  and  picked  for  detections  and  chirp  response  amplitudes.  The  chirp 
length  versus  epicentral  distance  data  points  so  obtained  are  plotted  in 
Figures  III-l,  III-2,  III-3,  and  III-4,  as  solid  circles.  In  these  figures, 
the  data  obtained  from  station -i  vents  detected  on  the  bandpass  filter  are 
denoted  by  open  circles.  TTie  least-inean-square-error  fit  shown  was  made 
on  the  data  points  represented  by  open  circles. 

B.  CHIRP  MATCHED  FILTER  RESULTS 

The  points  to  be  considered  in  the  evaluation  of  chirp  matched 
filters  are  dB  SNNR  improvements,  detection  level  improvements,  and 
surface-wave  magnitudes  derived  from  the  filter  responses.  These  will  be 
discussed  in  terms  of  data  from  the  two  regions. 
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CHIRP  LENGTH  VERSUS  DISTANCE  - CENTRAL  ASIA  - LR- 


CHIRP  LENGTH  VERSUS  DTSTAN'C 


Station-Events  Detected  on  Bandpass 


CHIRP  LENGTH  VERSUS  DISTANCE  - GREECE  - TURKEY 


1. 


dB  SNNR  Improvements 


The  dB  SNNR  improvement  of  the  chirp  matched  filter  over 
the  equivalent  (0.023-0.059  Hz)  bandpass  filter  was  calculated  by  the 
formula  given  in  Section  1-C.  SNNR  improvements  were  computed  for 
every  station-event  detected  on  the  bandpass  filter  response.  The  results 
are  listed  in  Table  III-2  for  the  central  Asia  region  and  Table  III-3  for  the 
Greece-Turkey  region. 

Considering  those  stations  which  detected  four  or  more  test 
events,  the  followir.g  comments  on  chirp  matched  filler  results  can  be  made. 
(The  term  "poor"  inaicates  a mean  SNNR  improvement  of  less  than  1 dB, 
"fair"  a mean  between  1 and  2 dB,  "good"  a mean  between  2 and  4 dB,  and 
"excellent"  a mean  greater  than  4 dB.) 

• Station  2 - SNNR  improvements  were  fair  for  test  events 
from  CENA. 

• Station  5 - SNNR  improvements  were  good  for  test  events 
from  CFNA  and  fair  for  test  events  from  GTUR. 

• Station  6 - SNNR  improvements  were  fair  for  test  events 
from  CENA  and  good  for  test  events  from  GTUR. 

• Station  8 - SNNR  i provements  were  good  for  test  events 
from  CENA. 

• Station  9 - SNNR  improvements  were  excellent  for  test 
events  from  GTUR. 

• Station  11  - SNNR  improvements  were  fair  for  test  ev  -its 
from  CENA. 
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TABLE  III-2 


CMF  SNNR  IMPROVEMENTS  FOR  CENTRAL  ASIA  EVENTS 

(LR-V) 

(PAGE  I OF  3) 


FVKNT  nn  SMNP  IMP.  nVFP 

VM1PER  STR'^TON  FOniVRLENT  PP  FILTER 
RR1  2 -O.U 

RRU  2 1.S 

RRS  2 1.S 

ROR  2 1.7 

qnq  2 2.0 

011  2 U. 7 ' 

om  2 1.1 

01S  2 1.1 

0 1 R 2 1 . R 

02R  2 1.2 

020  2 1.0 

OUR  2 0.7 

1001  2 U.R 

1002  2 2.2 

•1FAN  RNRP  If1PRnVERFNT=  1 . R2 
STRNnRpn  nEVTATjn*i=  l.io 


FVFNT 

DP  SNNR  IMP. 

OVER 

vniPFR 

station  fooival’=’nt  op 

FILTER 

RP1 

S 1.2 

PflU 

S 1.  R 

OIR 

S 1.7 

07fl 

S 3.R 

MEAN  SNNP 

IMPROVEMENT=  2.07 

OTA  NPAPn 

DEVTATTnN=  1.0P 
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CMF  SNNR  IMPROVEMENTS  FOR  CENTRAL  ASIA  EVENTS 


(LR-V) 

(PAGE  2 OF  3) 

FVFNT 

DR  SNNR  IMP.  nVFP 

ST  ATin*’ 

FOHTVALFNT  RP 

PR  1 

6 

2.0 

9F4 

6 

1.  3 

099 

6 

2.4 

9PR 

6 

0.9 

911 

f. 

1.6 

91f> 

fs 

4.9 

91P 

6 

2.7 

9?f 

6 

-0.9 

9?Q 

6 

1.9 

OSO 

6 

3.  1 

q92 

6 

2.9 

09U 

6 

2.7 

POO 

6 

-0.  1 

969 

6 

1.6 

070 

6 

1.4 

MEAN  SNNR 

T^PPOVEMENT=  2.03 

STANDAFP 

DFVT  »TT0N=  1.  3S 

FVFNT 

DP  SNNR  IMP. 

OVER 

NIIMPFP 

STATION 

ROOIVAIFNT  op 

RTLTRR 

001 

0 

1.6 

004 

0 

3.  4 

Q08 

0 

2.6 

P1 1 

0 

4.9 

910 

0 

6.9 

929 

0 

4.  1 

1 

046 

0 

2.  1 

990 

0 

9.  0 

96  9 

0 

-0.  4 

MS'IN  snnp 

T*’PnOV  FMFNT=  3.  «U 

STANDARD 

DEVIATION 

= 2.26 

III-ll 
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CMF  SNNR  IMPROVEMENTS  FOR  CENTRAL  ASIA  EVENTS 

<LR-V) 

(PAGE  3 OF  3) 


FVFNT 

PR  SNNR  TMR. 

OVER 

NH'inFR 

STATION 

FOniVAl.ENT  BP 

FIT  TEF 

qpu 

q 

ft.q 

POP 

q 

1.*^ 

1 

q 

q.u 

TIP 

q 

2.T 

qifl 

q 

6.  3 

qsp 

q 

2.S 

qsq 

q 

7.  R 

qp*^ 

q 

S.R 

q7R 

q 

ft.O 

•*EAN  SNNR 

T1PROVF'1FNT=  S.UI 

I standarp 

nEVTATTnN=  2.f>H 

FVFNT 

PB  SNNR  TIP.  OVER 

NIINpEP 

STAT'^ON 

roUIVATFNT’  PP  FIl.T’=’R 

qiR 

10 

4.3 

q2q 

10 

2.R 

BEAN  SWA'" 

TNPROVENFNT^  3.FB  | 

FVENT 

PR  SNNR  tmp,  over 

NUNRER 

ST’ATION 

EQHIVAtENT  RP  FILTER 

RP1 

1 1 

0.2 

RR4 

11 

1.B 

qm 

1 1 

3.P 

qip 

11 

0.6 

1 1 

O.R 

O^q 

1 1 

1.2 

qup 

1 1 

-0.6 

qfii 

1 1 

2.4 

qps 

1 1 

1.7 

nuu 

11 

-2.7 

lopq 

11 

2.P 

MEAN  SNNR  INPROVF?1ENT=  1.1^ 

STANPARP 

nFVTATTnN=  1.77 
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CMF  SNNR  IMPROVEMENTS  FOR  GREECE-TURKEY  EVENTS 

(LR-V) 

(PAGE  1 OF  2) 


FVRNT 

NnrRFP 

Rqn 

STATION 

2 

OP  SNNR  IMP. 
FQUT VALENT  RP 
U.2 

OVFP 

FTLTFR 

t-VpijT 

OR  SNNR  TIP. 

nVFP 

nii»inEP 

STATION 

FOniVALKNT  BP 

FITTEP 

R90 

S 

1.0 

911 

S 

1.0 

971 

S 

2.7 

9fll 

S 

2.  S 

10R1 

s 

o.u 

MFAN  SNNR  T'lPROVEIRNT^  1 . S2 
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li  is  immediately  obvious  that  more  data  is  needed  from  the 
Greece-Turkey  region  since  it  is  possible  to  compare  ihe  SNNR  improve- 
ments for  the  two  regions  at  only  three  stations.  At  station  5,  the  mean 
SNNR  improvement  for  CENA  is  higher  than  the  mean  SNNR  improvement 
for  GTUR.  At  station  6,  the  SNNR  improvements  for  GTUR  test  events 
have  a higher  mean  value  than  do  those  for  CENA  test  events.  At  station 
9i  the  SNNR  improvements  for  CENA  test  events  have  a higher  mean  value 
than  do  those  for  GTUR  test  events. 

The  standard  deviation  values  for  each  of  the  mean  SNNR 
improvements  indicate  there  is  considerable  variation  in  the  SNNR  improve- 
ments yielded  by  chirp  matched  filters.  For  example,  the  SNNR  improve- 
ment values  for  test  events  from  central  Asia  as  recorded  at  station  8 
range  from  -0.4  dB  to  6.  9 dB.  This  range  in  values  may  be  due  to  different 
source  mechanisms,  since  we  note  that  the  SNNR  improvements  can  be 
divided  roughly  into  two  groups,  one  of  low  SNNR  improvements  and  one  of 
high  SNNR  improvements.  Given  a sufficiently  large  data  base,  we  may  be 
able  to  determine  mean  SNNR  improvements  for  sub-regions  which  will 
have  low  standard  deviations. 

2.  Surface  Wave  Detection  Using  Chirp  Matched  Filters 

Using  the  detection  criteria  given  in  Section  II-C,  the 
detection  capability  of  chirp  matched  filters  was  evaluated  for  the  central 
Asia  region,  the  Greece- Turkey  region,  and  the  combined  central  Asia- 
Greece- Turkey  region.  This  evaluation  was  performed  in  terms  of  seismic 
events  and  the  VLPE  network.  The  results  are  as  follows: 

• CENA  - Fifteen  events  were  detected  on  the  bandpass  filter 

response.  An  additional  18  events  were  detected  on  the  chirp 
matched  filter  response,  resulting  in  a 120  percent  increase 
in  the  number  of  events  detected. 


GTUR  - Five  events  were  detected  on  the  bandpass  filter 
response.  An  additional  8 events  were  detected  on  the  chirp 
matched  filter  response,  resulting  in  a lf>0  percent  increase 
in  the  number  of  events  detected. 
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I 
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I 


f*' 


1 

I 

i 


Combined  CENA  and  GTUR  - Twenty  events  were  detected  on 
the  bandpass  filter  response.  An  additional  26  events  were 
detected  on  the  chirp  matched  filter  response,  resulting 
in  a 130  percent  increase  in  the  number  of  events  detected. 


Although  the  percentage  increase  in  the  number  of  events 
detected  differs  for  the  two  regions,  we  note  that  in  either  case  the  number 
of  detections  is  more  than  doubled  by  the  use  of  chirp  matched  filters. 

The  detection  capability  of  chirp  matched  filters  in  a network 

sense  is  illustrated  by  Figure  III-5,  which  shows  the  bodywave  magnitude 

distribution  and  maximum-likelihood  detectability  curve  for  the  combined 

central  Asia  and  Greece- Turkey  region.  Comparing  this  figure  with  the 

corresponding  bandpass  filter  maximum-likelihood  detectability  curve  of 

Figure  II-2,  we  find  that  the  use  of  chirp  matched  filters  lowered  the  50 

percent  detection  level  from  4.72  to  4.04  units  and  the  90  percent 

detection  level  from  5.  25  to  4.90  m units. 

b 

3.  Surface-Wave  Magnitude  From  Chirp  Matched  Filter  Data 

Surface-  wave  magnitudes  were  computed  from  chirp  matched 
filter  data  using  the  equation: 

“s  = f " I'OCio'DELTA)  +1.12 

This  is  the  same  as  the  equation  used  to  compute  M^  from  bandpass -filtered 

data  except  that  here  the  period  is  assumed  to  be  30  seconds,  since  no  period 

dB  . , . 

can  be  measured  from  a matched  filter  response.  The  term  — is  subtracted 
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to  remove  the  SNNR  improvement  due  to  the  matched  filter.  Since  the 
period  is  assumed  to  be  the  instrument  response  peak,  no  instrument  gain 
correction  is  needed  (the  instrument  response  is  flat  between  30  seconds 
and  40  seconds) . 

When  SNNR  improvement  values  were  available  for  four  or 
more  events  recorded  at  a station,  the  mean  of  those  values  was  used  in 
the  surface-wave  magnitude  calculation.  When  less  than  four  values  were 
available,  the  overall  average  SNNR  improvement  was  used. 

From  Tables  III-2  and  III-3,  the  SNNR  improvement  used 


SNNR  Improvements  j 

Statio  n 

CENA 

GTUR 

2 

1.8 

3.4* 

5 

2.  1 

1.5 

6 

2.  0 

3.3 

8 

3.4 

3.4* 

9 

5.4 

3.  0 

10 

2.  6* 

3.4* 

11 

1.  1 

3.  4* 

where  all  values  are  in  dB.  An  asterisk  denotes  the  use  of  the  overall 
average  SNNR  improvement  value. 

For  station-e''^ents  not  detected  on  the  bandpass  filter  response, 
the  sui  face-wave  magnitude  values  calculated  using  the  above  SNNR  improve- 
ments arc  shown  in  Figure  III-6.  When  M values  were  computed  at  two  or 

s 

more  stations  for  a given  event,  they  were  averaged  to  reduce  the  variance. 
The  resulting  plot  shows  network  average  values  with  one  or  more  stations 
reporting.  The  solid  line  is  a best  linear  fit  to  the  data,  computed  by 


III-18 


consiclc'r inj;  M and  to  be  independent  of  each  other  and  determining  a 

fit  by  minimizing  the  distances  normal  to  a line  and  the  data  points.  The 

method  of  computation  is  given  by  Lambert  et  al.  , 1973.  The  dashed  line 

is  a best  linear  fit  to  VLPE  bandpass  filtered  data  made  over  the*  range 

4.  2 < m ^5.5.  The  equation  ofthislineisM  =1.20m  -1.74 

b s b 

(Lambert  et  al.  , 197  3). 

The  four  circled  points  were  not  included  in  the  computation 

of  the  chirp  matched  filter  data  best  linear  fit  since  they  show  such  a large 

separation  from  the  rest  of  the  population.  Investigation  of  the  RMS  noise 

levels  for  the  days  on  which  these  ewents  occurred  indicate  that  these  high 

M values  are  due  to  high  noise  lc?vels.  For  c'xample,  consider  event  927 

from  the  GTUR  region.  This  event  had  an  m of  4.  6 and  an  M computed 

b s 

from  chirp  matched  filter  responses  of  3,88.  This  M value  is  an  average 

s 

of  th(‘  values  trieasured  at  stations  2 and  10.  At  station  2,  the  RMS  noise 

value  in  the  20-40  second  band  was  4,45  m/iand  the  computed  M value  was 

3.  30.  At  staticjn  10,  the  RMS  noise  value  in  the  same  band  was  28.  37m/i 

and  the-  computed  M was  4.46.  Thus,  it  appears  that  these  high  M values 

s s 

are  clue  at  lea' t in  part  to  high  RMS  noise  values  at  some  or  all  of  the 

stations  at  which  M measurements  were  made. 

s 

In  cjrder  for  an  event  which  was  not  detected  on  the  bandpass 
filter  response  tcj  be  detected  on  the  chirp  matched  filter  response,  the 
signal  cannot  lie  far  below  the  noise  level.  For  example,  assuming  a 6 dB 
SNNR  improvement  and  requiring  that  the  matched  filter  produce  a peak 
3 clB  above  the  noisc'  level  to  call  the  peak  a detection,  the  signal  cannot  be 
more  than  3 clB  below  the  noise  level  or  no  detection  will  occur.  With  this 
in  mind,  consider  Figure  III-6.  This  Figure  shows  that  as  the  m value 
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increases,  the  M values  of  events  detected  only  by  a matched  filter  depart 
s 

more  and  more  from  the  bandpass  filter  data  linear  fit  of  M 1.  20  m^  - 1.  74. 

At  high  m values,  we  are  above  the  50  percent  bandpass  filter  detection  level, 
b 

Thus,  for  high  rn  events,  the  only  surface  waves  not  detected  on  the  band- 

pass  filter  response  must  have  abnormally  low  surface  wave  magnitudes. 

Detections  by  chirp  filters  of  high  m^  events  which  were  not  detected  on  the 

bandpass  filter  response  show  that  this  is  the  case.  (Figure  III-6).  At  low 

values,  we  are  far  below  the  50  percent  bandpass  filter  detection  level,  and  the 

M values  for  events  detected  only  by  the  chirp  matched  filter  can  be  expected 
s 

to  be  close  to  the  M values  determined  from  bandpass  filtered  data  with  com- 

s 

parable  m,  values.  (We  note  that  in  the  case  of  an  event  detected  on  both 
b 

bandpass  and  chirp  matched  filter,  the  M value  measured  from  the  chirp 

response  closely  agrees  with  the  M value  measured  from  the  bandpiss  filter 

response.)  Thus,  M measurements  made  for  low  m events  detected  only  by 

s b 

chirp  filter  may  be  included  with  bandpass  data  in  - m plots. 

The  variance  in  the  data  set  may  be  ascribed  to: 

• Limited  data 

• No  control  of  the  m,  parameter 

b 

• No  accurate  period  estimate  of  the  signal 

• Station  bias 

• Path  bias 


Source  effects. 
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SECTION  IV 


REFERENCE  WAVEFORM  MATCHED  FILTER  EVALUATION 

A.  DISCUSSION 

Use  of  the  matched  filter  approach  implies  that  we  know  the 
expected  waveform  in  order  to  search  effectively  for  that  waveform  in  a 
noisy  record.  In  Section  III,  we  considered  the  chirp  matched  filter  method 
of  obtaining  the  expected  waveform.  We  now  turn  to  the  reference  waveform 
matched  filter  method  to  obtain  this  expected  waveform.  In  order  to  obtain 
a suitable  expected  waveform  and  assure  that  propagation  effects  were 
properly  accounted  for,  we  chose  the  surface  waves  from  an  event  detected 
on  the  bandpass  filter  which  had  an  epicenter  in  the  area  of  interest. 

The  criteria  for  choosing  a station-event  to  be  a reference 
waveform  were;  good  SNNR,  shallow  focus  (less  than  60  km),  and  location 
close  to  the  majority  of  the  other  events  of  the  region.  The  length  of  the 
reference  waveforms  was  chosen  in  the  following  manner:  for  events  at 
large  epicentral  distances,  the  length  was  selected  to  include  multipath 
energy,  since  small  changes  in  event  epicenter  location  would  not  be  expected 
to  significantly  change  the  multipath  structure.  This  situation  is  reversed 
for  events  at  small  epicentral  distances;  for  such  events,  small  changes  in 
event  epicenter  location  could  significantly  change  the  multipath  structure. 
Therefore,  the  lengths  of  reference  waveforms  having  small  epicentral 
distances  were  chosen  so  as  to  exclude  any  possible  multipath  energy. 

There  is  always  the  possibility  that  a given  refert.._e  wave- 
form will  perform  poorly  with  most  of  the  station-events  with  which  it  is 
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matched.  To  avoid  this,  two  reference  waveforms  were  chosen  for  each 
station- area  whenever  possible.  Both  reference  waveforms  were  then  matched 
with  all  api)ropriate  station-events.  The  better  of  the  two  was  tlu'n  selected 
primarily  on  the  basis  of  the  number  of  detections  and  secondarily  on  the 
average  SNNR  improvement  due  to  each,  h'or  this  reason,  the  variation  of 
SNNR  gain  due  to  the  RWMF  parameters  was  not  determined.  A list  of  the 
reference  waveforms  so  selected  s given  in  Table  IV- 1. 

B.  REFERENCE  WAVEEORM  MATCHED  FILTER  RESULTS 

As  in  the  preceding  section,  the  points  to  be  considc'red  in 
the  evaluation  of  reference  waveform  matched  filters  are  dB  SNNR  improve- 
ments, detection  level  improvements,  and  surface-wave  magnitudes  derived 
from  the  filter  responses.  These  points  will  be  discussed  in  terms  of  data 
from  the  two  regions  of  interest. 

1.  dB  SNNR  Improv<-ment 

The  dB  SNNR  improvement  of  the  reference  waveform 
matched  filter  over  the  equivalent  bandpass  filler  was  calculated  by  the 
formula  given  in  Section  II-C.  SNNR  improvement  value's  wc're  computed 
for  every  station-event  detected  on  the  bandpass  filter  response.  The 
results  are  listed  in  Table  IV-Z  for  the  central  Asia  region  and  Table  IV- 3 
for  the  Greece- Turkey  region. 

Considering  those  stations  which  cletec  ted  four  or  more*  test 
(■vents,  the  following  comments  on  reference  waveform  matched  filter 
results  can  be  made.  (The  terms,  "poor”,  "fair",  "good",  and  "(?xcellent" 
are  as  defined  in  Section  III-B-1.) 
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TABLE  lV-1 


LIST  OF  REFERENCE  WAVEFORM  MATCHED  FILTERS 


REGION 

1 

STATION 

RWMF 

CENA 

2 

881 

CENA 

5 

884 

CENA 

6 

926 

CENA 

8 

884 

CENA 

9 

884 

CENA 

10 

929 

CENA 

11 

884 

GTUR 

2 

890 

GTUR 

5 

983 

GTUR 

6 

890 

GTUR 

8 

890 

GTUR 

9 

897 

GTUR 

10 

897 

GTUR 

11 

890 
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• Station  2 - SNNR  improvements  were  poor  for  events  from 
central  Asia 

• Station  5 - SNNR  improvements  were  fair  for  events  from 
Greece-  Turkey 

• Station  6 - SNNR  improvements  were  fair  for  events  from 
central  Asia,  and  good  for  events  from  Greece- Turkey 

• Station  9 - SNNR  improvements  were  excellent  for  events 
from  central  Asia 

• Station  11  - SNNR  improvements  were  poor  for  events  from 
central  Asia. 

The  overall  average  SNNR  improvement  was  1 . 3 dB  for  central 
Asia,  and  3.  1 dB  for  Greece-Turkey.  A larger  data  base  for  Greece- Turkey 
will  be  necessary  to  allow  comparison  of  individual  stations  between  these 
regions  and  central  Asia. 

As  in  the  case  of  chirp  matched  filters,  the  standard  devia- 
tions associated  with  the  mean  values  of  SNNR  improvements  were  large. 

The  possibility  of  the  high  and  low  values  being  geographically  separated 
will  need  to  be  investigated  in  the  future. 

The  effect  of  reference  waveform-test  event  separation  upon 
SNNR  improvement  is  illustrated  bv  Figure  IV- 1.  (Only  the  CENA  data 
were  plotted,  since  there  were  too  few  GTUR  data  pc'ints  available  to  make 
a meaningful  plot.)  A straight-line  least-mean-square-error  fit  to  the 
data  points  is  shown. 
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Although  there  is  an  extremely  large  variation  in  the  data,  the  plot  seems  to 
indicate  that  the  SNNR  improvement  obtained  b/  the  reference  waveform 
matched  filters  decreases  gradually  with  increasing  reference  waveform- 
test  event  separation.  The  slope  of  the  fitted  line  is  -0.6  dB/1000  km,  as 
compared  to  a slope  of  -3.0  dB/1000  km  for  ALFA  data. 

2.  Surface-Wave  Detection  Using  Reference  Waveform  Matched 

Filters 

Using  the  detection  criteria  given  in  Section  II-C,  the 
detection  capability  of  reference  wnveform  matched  filters  was  evaluated 
for  the  central  Asia  region,  the  Greece- Turkey  region,  and  the  combined 
central  Asia-Greece  - Turkey  region.  This  evaluation  was  performed  in 
terms  ol  seismic  events  and  the  VLPE  network.  The  results  are  as  follows: 

• CENA  - Fifteen  events  were  detected  on  the  bandnass  filter 
response.  An  additional  20  events  were  detected  on  the 
reference  waveform  matched  filter  response,  resulting  in 

a 133  percent  increase  in  the  number  of  events  detected. 

• GTUR  - Five  events  were  detected  on  the  bandpass  filter 
response.  An  additional  8 events  were  detected  on  the 
reference  waveform  matched  filter  response,  resulting  in 
a 160  percent  increase  in  the  number  of  events  detected. 

• Combined  CENA  and  GTUR  - Twenty  events  were  detected  on 
the  bandpass  filter  response.  An  additional  28  events  were 
resulting  in  a 140  percent  increase  in  the  number  of  events 
detected. 
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Although  the  percentage  increase  in  the  number  of  events 
detected  differs  for  the  two  individual  regions,  we  note  that  overall  the 
number  of  detections  is  more  than  doubled  by  the  use  of  reference  wave- 
form matched  filters. 

The  detection  capability  of  reference  waveform  matched 

filters  in  a network  sense  is  illustrateo  by  Figure  IV-2,  which  shows  the 

bodywave  magnitude  distribution  ana  maximum-likelihood  detectability 

curve  for  the  combined  central  Asia  and  Greece- Turkey  region  events. 

Comparing  this  figure  with  the  corresponding  bandpass  filter  maximum- 

likelihood  detectability  curve  of  Figure  II-2,  we  find  that  the  use  of 

reference  waveform  matched  filters  lowered  the  50  percent  detection  level 

from  4.  72  to  4.01  m units  and  the  90  percent  detection  level  from  5.  25 

b 

to  4.91  m,  units, 
b 

3.  Surface-Wave  Magnitudes  from  Reference  Waveform  Matched 

Filter  Data 

Surface-wave  magnitudes  were  computed  from  reference 
waveform  matched  filter  data  using  the  method  described  in  Section  II-B-3 
for  those  station -events  detected  on  the  matched  filter  response  but  not  on  the 
bandpass  filter  response.  When  SNNR  improvement  values  were  available 
for  four  or  more  events  recorded  at  a station,  the  average  of  those  values 
was  used  in  the  surface-wave  magnitude  calculation.  When  less  than  four 
values  were  available,  the  overall  average  SNNR  improvement  was  used. 

From  Tables  IV-2  and  IV-3,  the  SNNR  improvements  used 


were: 
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FIGURE  IV -2 

DETECTION  STATIST  ICS  FOR  THE  COMBINED  REGION-RWM} 
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where  all  values  are  in  clB.  An  asterisk  denotes  the  use  of  the  overall 
average  SNNR  improvement  value. 
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The  surface-wave  magnitude  values  calculated  using  the 
above  SNNR  improvement  values  are  shown  in  Figure  IV- 3.  When 
valuf's  were  computed  at  two  or  more  stations  for  a given  event,  they  were 
averaged  to  reduce  the  variance.  The  resulting  plot  shows  network  average 

M values  with  one  or  more  stations  reporting.  Tin*  M values  are  for 
s ^ 

events  detected  by  reference  waveform  matched  filters  but  not  by  the  band- 
pass filtered  data  as  described  in  Section  III-B-3.  The  solid  line  is  a best 
linear  fit  to  the  data,  computed  by  considering  M and  m to  be  independent 
of  each  other  and  determining  a fit  by  minimizing  the  distances  normal  to  a 
line  and  the  data  points.  The  circled  points  were  not  included  in  the  computa- 
tion of  this  linear  fit,  since  three  show  anomalously  high  values  and  one  an 
anomalously  low  value.  As  in  the  case  of  the  chirp  matched  filler  data,  it 
is  believed  that  the  high  values  are  due  at  least  in  part  to  high  RMS  noise 

values.  The  anomalously  low  M value  may  be  a false  alarm. 

s 

The  data  of  Figure  IV- 3 present  the  same  picture  as  did  the 
chirp  matched  filter  data.  At  low  m values,  the  M values  computed  for 
events  detected  only  by  the  reference  waveform  matched  filter  s arc  com- 
parable to  the  M values  for  bandpass  filtered  data.  At  higher  m values 
s " 

(near  or  above  the  50  percent  detection  level),  the  values  computed  for 

events  detected  only  by  the  reference  waveform  matched  filters  are  much 
lower  than  the  M values  determined  from  bandpass  filtered  data. 
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SECTION  V 


THREE-COMPONENT  ADAPTIVE  PROCESSOR 
EVALUATION 


It  has  been  previously  reported  (Lane,  1973)  that  the  three- 
component  adaptive  processor  can  yield  as  much  as  8 dB  improvement  in 
the  signal -plus-noise  to  noise  ratio  when  applied  to  single  site  data.  However, 
due  to  unmatched  instrumental  phase  responses  between  the  horizontal  and 
vertical  components  of  the  VLPE  stations,  it  is  not  expected  that  such  SNNR 
improvements  will  be  realized  when  this  method  is  applied  to  VLPE  data. 

In  the  following,  we  assess  the  value  of  the  three-component  adaptive  pro- 
cessor when  applied  to  VLPE  data  in  its  present  form  to  determine  whether 
this  method  can  be  used  in  spite  of  this  problem. 

A.  DISCUSSION 

The  three-component  adaptive  filter  (TCA)  is  designed  to 
improve  the  detectability  of  long-period  Love  and  Rayleigh  waves  by  utilizing 
the  known  phase  relationships  between  these  surface  waves.  (Lane  1973). 

The  process  examines  three  mutually  perpendicular  seismometer  traces  and 
after  dividing  the  data  into  overlapping  segments  designs  a new  filter  for 
each  segment  (i.  e.  , it  is  adaptive). 

When  considering  Rayleigh  waves,  the  TCA  examines  only  the 

radial  and  vertical  components.  Ideally,  fundamental  Rayleigh  wave  signals 
o 

are  90  out  of  phase  on  these  components,  and  consequently,  the  filter  is 
designed  to  emphasize  this  apparent  polarization  of  particle  motion. 


V-1 


Consider  the  k-th  segment  of  data,  and  let  V'-  ( u } and  R’l'  ( u ) 

k k 

be  the  Fourier  transforms  at  the  l/-th  frequency  of  the  vertical  and  radial 
traces,  respectively.  As  complex  functions,  these  components  may  be 
expressed  as; 


iv)  = 

k 


V»  (V)  c 
k 


R^k)  <-■ 


iO  (K) 
k 


i‘/>^(l/  ) 


where  i = V^l  and  (j^(  V),  V ) are  the  phase  spectra  of  the  vertical  and 

radial  traces,  respectively.  Then,  the  relative  phase  angle  tj/  ( p)  is  given  by; 

Ic 

( l^)  = 0 ( I/)  - (/>  (p ). 

k k k 

and  the  filter  transfer  function  (filter  weight  lunction)  is  defined  by; 

U ) - sin^  (ijj^iv)). 


where  N must  be  an  even  integer  in  order  to  guarantee  that  F ( P ) is  non- 

k 

negative.  This  exponent  determines  the  degree  of  rejection  by  the  filter  for 
frequencies  whose  relative  phase  angles  are  undesirable.  Since  Rayleigh 
wave  signals  are  usually  contaminated  by  noise  (and,  occasionally  by  higher 
modes),  the  relative  phase  angle  will  never  be  exactly  90°,  and  subsequently, 
care  must  be  taken  to  not  use  too  large  an  exponent  or  else  significant  amounts 
of  signal  energy  will  be  rejected.  For  the  data  processed  herein,  a value  of 
N 6 was  selected  with  segments  ciiosen  to  be  64  points  (128  seconds)  in 
length.  These  values  were  suggested  by  Lane,  (1973).  Other  values  are 
currently  being  examined  to  determine  whether  or  not  these  are  the  optimum 
values. 

The  bandpassed  Fouri'^r  transforms  of  the  traces  are  then 

multiplied  by  this  filter  to  give  the  transforms  V ( l^),  R ( p ) of  the  filtered 

k V 


segments; 


I 


I*  i 
I I 


i i 


V (I/)  = ( I/)  • F f u) 

k k k 

R,(l^)  = R-  (l^)  - F ( U). 
k k k 

These  traces  are  then  inverse  transformed  back  to  the  time  domain.  Next, 
the  k + 1 segment  of  time  domain  data  is  examined  and  defined  to  begin  at 
the  mid-point  of  the  k-th  segment.  After  filtering  as  above,  a weighted 
average  of  the  resultant  k and  k + 1 segments  is  formed  over  their  common 
interval.  In  particular,  at  the  i-th  time  point  the  final  value  is  given  by; 

\<V-  Vi'V<f>+ W'-y’' 

wher  e 

/ = length  of  overlapping  segment, 

d.  = distance  from  beginning  of  k + 1 segment  to  t -th  point. 

1 i 

When  considering  Love  waves,  the  TCA  analyses  the  transvc'rse 

and  radial  components  and  enhances  those  signal  elements  which  are  on 

azimuth.  In  particular,  the  azimuthal  filter  is  defined  (similarly  to  the 

N 

polarization  filter)  by;  F ( u ) = cos  ( ^ ( l^)).  where  N is  an  even  integer, 

K K 

U is  the  frequency  and  ^ ( U ) is  the  difference  between  the  expected  azimuth 
of  arrival  of  the  suspected  event,  o , and  the  apparent  azimuth  of  arrival  as 
determined  by  the  transverse  and  radial  traces  for  the  k-th  segment.  In 
particu’ar , 


i , ( P)  = o - tan 
k 


-1 


where  R’:'  (p;and  T-!'  ( p ) are  the  Fourier  transforms  at  the  P-th  frequency 
for  the  k -th  segment  of  the  radial  and  the  transverse  traces,  respectively. 

The  transverse  and  radial  time  domain  traces  are  then  treated 
analogously  to  the  Rayleigh  wave  process  presented  above  but  with  the  trans- 
verse component  being  considered  rather  than  the  vertical  component  and 
with  the  filter  transfer  function  redesigned  as  above. 
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The  data  base  used  to  evaluate  the  TCA  technique  was 
essentially  the  same  as  that  used  to  evaluate  the  matched  filter  techniques. 
However,  if  one  component  of  a station-event  contained  malfunctions,  it 
could  not  be  processed  by  the  TCA,  but  could  be  processed  by  matched 
filters.  Therefore,  the  data  base  for  the  TCA  evaluation  is  somewhat 
smaller  than  that  used  for  the  matched  filter  evaluation.  From  the  Tables 
of  the  Appendix  A,  we  see  that  204  station- events  from  central  Asia  and  97 
station- events  from  Greece- Turkey  were  available  for  processing  by  the 
TCA  technique. 

b.  THREE-COMPONENT  ADAPTIVE  PROCESSOR  RESULTS 

As  in  tlie  preceding  sections,  the  points  to  be  considcrco  in 
the  evaluation  of  the  TCA  technique  are  dB  SNNR  improvements,  detection 
level  improvements,  and  surface-wave  magnitudes  derived  from  the  TCA 
responses.  These  points  will  be  discussed  in  terms  of  data  from  the  two 
regions  of  interest. 

1.  dD  SNNR  Improvement 

The  dB  SNNR  improvement  of  the  TCA  processor  over  the 
equivalent  bandpass  filter  was  calculated  by  the  formula  given  in  Section  II-C. 
The  RMS  noise  values  v'cre  determined  from  noise  gates  immediately  pre- 
ceding the  signal  gates.  The  results  are  listed  in  Table  V-1  for  central  Asia 
and  Table  V-2  for  Greece- Turkey . 

Considering  those  stations  which  detected  four  or  more  test 
events,  the  following  comments  on  the  three-component  adaptive  processor 
results  can  be  made.  (The  terms  "poor",  "fair",  "good",  and  "excellent" 
are  defined  in  Section  III-B-l.) 
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• Station  2 - SNNR  improvements  were  fair  for  events  from 
CENA 

• Station  6 - SNNR  improvements  were  poor  for  events  from 
CENA  and  GTUR 

• Station  8 - SNNR  improvements  were  fair  for  events  from 
CENA 

• Station  9 - SNNR  improvements  were  fair  for  events  from 
CENA  and  poor  for  events  from  GTUR 

• Station  11  - SNNR  improvements  were  good  for  events  from 
CENA. 

The  overall  average  SNNR  improvement  was  1.6  dB  for  central 
Asia  and  0.  6 dB  for  Greece- Turkey.  Although  a larger  data  base  for  Greece- 
Turkey  would  be  necessary  to  allow  comparison  of  more  individual  stations, 
it  is  immediately  obvious  that  only  low  mean  values  of  SNNR  improvement 
were  obtained  from  this  technique.  These  low  values  of  SNNR  improvement 
are  primarily  due  to  unmatcned  instrumental  phase  res})onses  between  the 
horizontal  and  vertical  components.  Lambert  (Lambert  et  al.  , 1973)  recog- 
nized this  problem  and  discussed  it  in  terms  of  the  errors  in  analytic  rotation 
of  the  north  and  east  components  to  transverse  and  radial  components.  Visual 
observation  of  bandpass  filtered  data  shows  that,  in  general,  the  vertical  and 
horizontal  Rayleigh  waves  are  not  90  out  of  phase. 

2.  Surface  Wave  Detection  Using  the  Three-Component  Adaptive 

Processor 

Using  the  detection  criteria  of  Section  II-C,  the  detection 
capability  of  the  three-component  adaptive  processor  was  evaluated  for  the 
central  Asia  region,  the  Greece- Turkey  region,  and  the  combined  central 
Asia  and  Greece  - Turkey  regions.  This  evaluation  was  performed  in  terms  of 
seismic  events  and  the  VLPE  network.  The  results  are  as  follows: 
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• CENA  - Fifteen  events  were  detected  on  the  bandpass  filter 
response.  An  additional  two  events  were  detected  on  the 
three-component  adaptive  processor  response,  resulting 

in  a 1 3 percent  increase  in  the  number  of  events  detected. 

• GTUR  - Five  events  were  detected  on  the  bandpass  filter 
response.  No  additional  events  were  detected  on  the  three- 
component  adaptive  processor  response. 

• Combined  CENA  and  GTUR  - Twenty  events  were  detected 
on  the  bandpass  filter  response.  An  additional  two  events 
were  detected  on  the  three-component  adaptive  processor 
response,  resulting  in  a 10  percent  increase  in  the  number 
of  -vents  detected. 

Thus,  we  see  that  the  three-component  adaptive  processor, 
as  currently  used,  does  not  make  any  significant  changes  in  the  number  of 
events  detected. 

The  detection  capability  of  the  three-component  adaptive  pro- 
f’essor  is  illustrated  by  Figure  V-1.  This  figure  shows  the  bodywave  mag- 
nitude distribution  and  maximum-likelihood  detectability  curves  for  the 
combined  central  Asia  and  Greece  - Turkey  region  events.  Comparing  this 
figure  with  the  corresponding  bandpass  filter  maximum-likelihood  detecta- 
bility curve  of  Figure  II-2,  we  find  that  the  50  and  90  percent  detection  levels 
are  virtually  unchanged. 

3.  M Computations  for  TCA 

s 

Within  the  TCA  program  an  M value  for  the  vertical  component 
was  computed  for  the  maximum  peak-to-peak  amplitude  within  a signal  gate. 
When  the  event  was  considered  to  be  detected,  this  was  its  assigned  (initial) 

M value,  designated  M . The  following  program  formula  was  employed; 
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where: 


A 

T 


(t.) 


) + log  A + 1.12 


T * Y (t  ) 
k 1 


= Maximum  peak-to-peak  amplitude  (m/i) 

= Period  corresponding  to  A 

= Epicentral  distance  from  event  to  station  (degrees) 

= Final  filter  weight  utilized  at  the  t -th  time  point  where 

i 

the  t. -th  point  was  the  mid-zero  crossing  of  the  measured 
peak-to-peak  wave. 


This  formula  was  subsequently  corrected  for  the  station's  instrument 

response  to  arrive  at  a final  M value,  M (TCA): 

s s 

M (TCA)  - m'  - log  (G  (T)) 
s s 10 

where  G ( T)  = station's  instrument  response  at  period  T. 

In  an  effort  to  assure  that  reasonable  M values  were  being 

s 

generated  by  the  program  we  examined  certain  February  197  3 events.  These 
events  were  selected  to  fulfill  the  following  requirements; 

• The  events  were  detected  by  the  bandpass  filter  on  both 
vertical  and  transverse  components  (at  each  station  considered.  ) 

• All  three  components  were  functioning.  (I’his  insured  that 
vertical  and  transverse  TCA  output  would  be  meaningful). 

• At  least  two  stations  detected  the  event  and  these  stations 

were  on  approxim.ately  the  same  epicenter-to-station  azimuth. 

This  requirement  minimized  event  M (bandpass)  variations 

s 

resulting  from  path  effects. 
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There  were  six  station -events  which  satisfied  all  of  the  above 

requirements.  For  these  events  the  M (bandpass)  and  the  M (TCA)  were 

s s 

found  to  be  comparable.  The  results  are  presented  in  Table  V-3.  Thus, 

we  concluded  that  the  M (TCA)  values  produced  by  the  program  were 

s 

reasonable. 

As  shown  in  Figure  V-2,  the  poor  detection  capability  of  the 

three-component  adaptive  processor  resulted  in  only  12  M values  of  events 

detected  on  the  three-component  adaptive  processor  but  not  the  bandpass  filter. 

For  comparison,  the  linear  fit  to  VLPE  bandpass  filtered  data  (Lambert 

et  al.  , 197  1)  is  presc-ntc'd.  Since  th(*re  arc'  only  12  M values  available  in  a 

s 

narro  v 'll  range,  a linear  fit  was  not  computed  for  these  values.  There  is  no 
b 

trc'nd  obsc'rvable  in  thc'sc'  data  points  comparable'  to  that  of  the  matched  filter 

data  of  Figure  III-6  and  IV-1.  It  is  not  known  whether  this  is  due  to  errors 

in  the  computation  of  M from  three-component  adaptive  processor  data  or 

s 

is  merc'ly  clue  to  the  limiti-d  number  of  data  points. 
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TABLE  V-3 

M (TCA)  VERSUS  M (BANDPASS) 
s s 


Station 


4.  32 


M (TCA)  (Bandpass) 
s 


3.98 
4.  20 


4.04 
3.  25 


3.  82 


SECTION  VI 


DISCUSSION  AND  SUMMARY 

A.  COMPARISON  OF  THE  TECHNIQUES 

Each  of  the  three  data  enhancement  techniques  has  now  been 
evaluated  in  terms  of  dB  SNNR  improvement,  detection  level  improvc'ment, 
and  surface-wave  magnitudes.  We  now  compare  these  techniques  in  an 
attempt  to  determine  wliich,  if  any,  is  superior  to  the  others  for  a given  region. 

1.  Comparison  of  SNNR  Improvements 

The  mean  SNNR  improvements  and  associated  standard 
deviations  for  each  technique  arc  listed  by  region  in  Table  VI-1.  In  each 
region,  the  best  technique  is  judged  to  be  that  one  which  displays  the 
largest  mean  SNNR  improvement  with  the  smallest  associated  deviation. 

Using  this  criterion,  the  best  technique  for  data  enhancement,  by  region  is: 

CENA  - The  chirp  matched  filter  technique  outperformed  the  other 
two . 

GTUR  - The  chirp  matched  filter  technique  outperformed  the  other 
two.  However,  the  reference  waveform  matched  filter 
technique  was  almost  as  good. 

Another  comparison  of  the  three  techniques  is  given  by 
Figures  VI-1  and  VI-2.  Figure  VI-1  indicates  that  approximately  75  percent 
of  the  Central  Asia  test  events  showed  higher  chirp  SNNR  improvements 
than  reference  waveform  SNNR  improvements,  and  approximately  67  percent 


TABLE  VI-1 

COMPARISON  OF  CMF,  RWMF,  AND  TCA  SNNR  IMPROVEMENTS 


Region 

Technique 

Mean 

Standard 

Deviation 

CENA 

CMF 

2.  6 

2.  2 

RWMF 

1.  3 

2.9 

TCA 

1.  6 

2.  3 

GTUR 

CMF 

3.4 

2.  0 

RWMF 

2.  8 

2.  2 

TCA 

0.  6 

1.  9 
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FIGURE  VI-l 

CUMULATIVE  DISTRIBUTION  OF  CMF-RWMF  SNNR  IMP  (A) 
AND  CMF-TCA  SNNR  IMP  (B)  FOR  CENA  TEST  EVENTS 
(VERTICAL  COMPONENT) 
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FIGURE  VI-2 

CUMULATIVE  DISTRIBUTION  OF  CMF-RWMF  SNNR  IMP  (A) 
AND  CMF-TCA  SNNR  IMP  (B)  FOR  GTUR  TEST  EVENTS 
(VERTICAL  COMPONENT) 
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Percentage  of  Events  Percentage  of  Events 


of  the  central  Asia  test  events  showed  higher  cliirp  SNNR  improvements 
than  threc-component  adaptive  processor  SNNR  improvements.  Figure  VI-2 
indicates  that  Greece- Turkey  test  events  yielded  about  tlie  same  SNNR 
improvc-ments  on  chirp  and  reference  waveform  matclied  filters,  while 
approximately  85  percent  of  the  test  events  from  tliis  region  showrul  higher 
chirp  SNNR  improvements  than  three -component  adaptive  processor  SNNR 
improvements.  Th(>refore,  by  this  criterion,  chirp  matched  filters  out- 
performed both  the  reference  waveform  and  three- -component  adaptive  pro- 
cessor on  central  Asia  events.  For  Greece- Turkey  events,  chir]>  and 
reference  waveform  matched  filters  out-performed  the  three-component 
adaptive  filter  by  tlie  same  amount.  This  is  in  agreement  with  the  preceding 
judgment  made  on  the  basis  of  mean  SNNR  improvement  and  associated 
standard  deviation. 

2.  Gomparison  of  Detection  Capabilities 

The  second  point  of  comparison  is  the  relative  ability  of  the 
three  tccliniques  to  detect  signals  which  w'ere  not  detected  on  the  bandpass 
filt(‘r  response.  For  tlie  comljined  region,  use  of  chirp  matched  filters 
resulted  in  a 130  percent  increase  in  the  number  of  events  detected  and  use 
of  reference  waveform  matched  filters  resulted  in  a HO  percent  increase 
in  the  number  of  events  detected.  The  use  of  the  three -component  adaptive- 
processor  resulted  in  only  a 10  percent  increase  in  the  number  of  events 
detected.  Thus,  in  terms  of  tlie  increase  in  the  number  of  events  detected, 
use  of  chirp  or  reference  waveform  matched  filters  mort-  than  doubled  the 
number  of  detections,  while  use  of  the  three-component  adaptive  processor 
resulted  in  very  little  improvement  in  the  number  of  detections. 

An  increase  in  the  number  of  detections  results  in  a cU-creasi- 
in  the  maximum4ikelihood  detection  levels.  This  is  illustrated  by  I’ahle  VI-2, 
which  lists  the  50  and  90  ptircent  detection  levels  resulting  from  application 
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TABLE  VI-2 

DETECTION  LEVELS  DUE  TO  USE  OF  THE  DATA 
ENHANCEMENT  TECHNIQUES 


Technique 


m Detection  Levels 
b 


50% 


90% 


Bandpass  Filter 

CMF 

RWMF 

TCA 


4.  72  + 0.  09 
4.  04  + 0.  11 
4.  01  + 0.  11 
4.68  + 0.  10 


5.  25  + 0.  17 

4.  90  + 0.  20 

4.91  + 0.  21 

5.  31  + 0.  10 


of  the  maximum-likelihood  method  to  each  of  the  three  data  enhancement 

te  hniquer  data  sets.  The  bandpass  filter  detection  levels  are  listed  as  a 

reference.  All  values  arc  for  the  combined  region.  This  table  indicates 

that  use  of  chirp  or  reference  waveform  matched  filters  decreases  the  50 

percent  detection  level  by  approximately  0.7  m units  and  the  90  percent 

detection  level  by  approximately  0.3  m units.  The  changes  in  the  50  and 

b 

90  percent  detection  levels  due  to  use  of  the  thrco-componcnt  adaptive 
filter  arc  too  sn-^all  to  be  significant. 

It  has  been  noted  previously  (Lane,  1973)  tliat  an  SNNR  improve- 
ment of  6 dB  implies  a reduction  of  about  0.  ' in  the  bodywave  magnitude  at 
wliich  50  percent  of  all  events  are  detected  (the  50  percent  detection  level). 

This  in  turn  implies  a floubling  of  the  total  number  of  events  detected.  In 
this  report,  we  have  noted  that  use  of  cliirp  or  reference  waveform  matclied 
filters  results  in  SNNR  improvements  of  about  3 dB,  a reduction  of  about 
0.7  in  the  50  percent  detection  level,  and  a factor  of  about  '*.4  increase  in 
the  total  number  of  events  detected.  The  explanation  for  these  apparent 
anomalies  is  as  follows.  First,  tlie  low  mean  SNNR  improvements  (It'ss 
than  2 dB)  computed  at  some  stations  probably  do  not  represent  the  SNNR 
improvement  produced  when  a station-event  which  was  not  detected  on  tlie 
bandpass  response  becomes  visible  on  a matched  filter  response.  This  must 
be  true  for  central  Asia  station -events  of  this  type  whicli  are  detected  by 
reference  waveform  matched  filters  at  Stations  2 and  11,  since  the  mean 
SNNR  improvements  in  these  cases  are  negative,  T)ius,  t)ie  SNNR  improve- 
ments are  probably  higher  than  tlie  approximate  3 dB  improvement  previously 
computed.  Next,  we  recall  that  an  event  is  considered  to  be  delected  only 
if  it  is  detected  at  two  or  more  stations.  Those  events  wbicli  were  detected 
on  the  bandpass  filter  response  at  only  one  station  (a  total  of  19  events)  were 
therefore  listed  as  non-detected  events.  Therefore,  a detection  by  a matclied 
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filter  at  only  one  station  other  than  the  one  at  which  it  was  detected  on  band- 
pass will  change  the  detection  status  from  non-detected  to  detected.  Since 
it  is  more  likely  to  detect  an  event  using  matched  filters  at  one  station  than 
at  two  for  a given  SNNR  improvement,  we  see  that,  for  the  VLPE  network, 
it  is  possible  to  more  than  double  the  number  of  detections  using  matched 
filters  when  the  mean  SNNR  improvement  is  less  than  6 dB. 


The  factor  of  2.  4 increase  in  the  total  number  of  events 

detected  due  to  use  of  chirp  or  reference  waveform  matched  filters  implies 

a reduction  of  0.4  m units  in  the  50  percent  detection  level.  For  this 

b 

particular  data  base,  however,  we  have  a 0.  7 m unit  reduction  m the  50 

b 

percent  detection  level.  To  resolve  this  anomaly,  we  must  have  a larger 
data  base,  especially  at  the  lower  bodywave  magnitude  values. 


Since  the  chirp  and  reference  waveform  matched  filters  yield 
approximately  the  same  improvement  in  detection  levels,  let  us  consider  what 
happens  if  the  detection  results  are  combined.  In  addition  to  the  detection 
criteria  used  previously,  an  event  is  now  considered  to  be  detected  if  it  is 
detected  on  either  the  chirp  or  reference  waveform  matched  filter  resf>onse. 
Twenty  of  the  79  events  of  the  combined  region  were  detected  on  the  bandpass 
filter  response.  Use  of  this  new  detection  criterion  results  in  a total  of  35 
additional  detections.  The  result  of  this  is  shown  in  Figure  VI-3,  which 
indicates  that  under  this  criterion,  the  50  percent  detection  level  is  3.74  + 
0.12  and  the  90  percent  detection  level  is  4.51  ^ 0.17  for  the  combined  region. 
Thus,  the  detection  levels  are  lowered  significantly  (about  0.  3 to  0.  4 m 

b 

units)  relative  to  those  where  detection  is  by  chirp  or  reference  waveform 
matched  filter  alone. 


i 


( 
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FIGURE  Vl-3 

DETECTION  STATISTICS  FOR  THE  COMBINED  REGION  REQUIRING 

DETECTION  ON  C MF  OR  R WMF 
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3. 


Comparison  of  Surface-Wave  Magnitudes 


The  point  of  comparison  for  surface-wave  magnitudes  derived 

from  each  of  the  techniques  is  the  linear  fit  made  to  each  set  of  M : m data, 

s ■ b 

where  the  M values  are  for  events  not  detected  on  the  bandpass  filter.  For 
s 

the  combined  region,  the  equations  for  these  linear  fits  are: 


M = 

0.  25  m 

i 

1.  94 

for 

CMF  data 

s 

b 

M = 

0.  28  m 

i 

1.  81 

for 

RWMF  data 

s b 

(No  linear  fit  was  computed  for  the  TCA  data,  since  only  12  data  points  were 

available.)  From  these  equations,  we  sec  that  the  two  matched  filter  methods 

produce  values  for  surface-wave  magnitudes  which  have  essentially  the  same 

M : m,  relationship.  The  small  differences  between  the  two  relationships 
s b 

are  mostly  due  to  inaccuracies  in  the  dD  SNNR  improvements. 

To  compare  M values  computed  from  TCA  data  with  those 
s 

computed  from  matched  filler  data,  we  consicl-  r Table  VI- 3,  which  is  a con- 
densed version  of  Table  B-1  of  Appendix  B.  Table  VI-3  lists  only  those 

events  for  which  an  M value  was  computed  from  the  CENA  three-component 

s 

adaptive  processor  data.  This  table  shows  that  all  but  one  of  the  TCA  surface - 

wave  magnitudes  are  higher  than  cither  of  the  corresponding  matched  filter 

surface-wave  magnitudes.  Since  it  is  beli>  ved  that  the  matched  filter  M 

values  are  representative  of  the  events  detected,  it  will  be  necessary  in 

future  work  to  re-examine  the  manner  in  which  M is  computed  from  TCA  data. 

s 

B.  SUMMARY  OF  RESULTS 

The  major  conclusions  of  this  evaluation  of  the  chirp  matched 
filter,  rc'fercnce  waveform  matched  filter,  and  three -component  adaptive 
processor  data  enhancement  techniques  are: 
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TABLE  VI-3 

COMPARISON  OF  SURFACE-WAVE  MAGNITUDES 


Event 

Number 

""b 

Surface-Wave  Magnitude 

CMF 

RWMF 

TCA 

885 

4.  8 

3.  13 

2.  69 

3.  03 

887 

4.  7 

2.  72 

2.  60 

3.  07 

914 

4.  6 

3.03 

3.  39 

4.17 

955 

4.  2 

3.  37 

3.  38 

4.  47 

961 

4.  3 

3.08 

3.04 

3.  52 

1053 

5.  0 

3.  10 

- 

3.  69 

1084 

4.  5 

- 

3.03 

4.  22 

« 


In  the  two  seismic  regions,  the  chirp  matched  filter 
technique  outperformed  the  other  two  techniques  in  terms 
of  overall  mean  LR  SNNR  improvement.  (Since  the  standard 
deviations  of  the  mean  improvements  were  large,  it  is  not 
meaningful  to  attempt  a quantitative  statement  of  relative 
performance. ) 


• Even  though  the  overall  mean  LR  SNNR  improvement  for  a 

given  technique  applied  to  events  from  a given  region  may 
be  low,  the  improvement  in  detection  may  be  good. 

® In  terms  of  the  increase  in  the  number  of  events  detected, 

the  two  matched  filter  techniques  performed  equally  well 
and  far  outperformed  the  three-component  adaptive  processor 
technique. 


In  terms  of  the  detection  level  improvement  of  the  network 
considered  in  this  report,  the  two  matched  filter  techniques 
performed  equally  well  and  far  outperformed  the  three-com- 
ponent adaptive  processor  technique.  When  applied  to  the 

data  set  of  this  report,  both  yielded  a 0.  7 m,  unit  reduction 

b 

in  the  50  percent  twe -station  detection  level  and  a 0.  3 m 

b 

unit  reduction  in  the  90  percent  detection  level. 


When  dealing  with  M values  of  events  detected  only  by  a 

matched  filter,  values  comparable  to  those  from  bandpass 

filtered  data  can  be  expected  for  m,  values  below  the  50  per- 

b 

cent  bandpass  filter  detection  level.  For  higher  m values,  the 

b 

M values  for  events  detected  only  by  a matched  filter  can  be 
s 

expected  to  be  much  lower  than  ’.he  M values  determined  from 


s 

bandpass  filtered  data  with  comparable  m,  values. 

b 
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• Ovc-rall,  thcrt'  is  no  clear  superiority  of  one  matt  heel  filter 
technique  over  the  other  for  the  se  of  stations  considered 
in  this  report.  Both  are  superior  to  the  three-component 
adaptive'  jirocessor  techniejue  as  it  is  presently  use'el. 

• The  poor  performance  of  tlie  three-component  adaptive  pro- 
cessor is  ehie  not  tt)  some  intrinsic  flaw  in  the  methoel  but 

to  the'  unmatched  instrumental  phase  re'sponse's  betwe'en  the 
horizontal  anel  vertical  compements  of  tlie  VLPE  stations. 

C.  FUTURE  WORK 

The  followinj;  points  shoulel  be'  consielere-el  in  any  future'  work 
along  the  lines  of  this  report: 

• The  elata  base*  should  be  increase'el  --  i-*iore  e've'nts  anei  stations 
should  be  inve'stigateel  to  be-tter  assess  tlie  capabilitie-s  of 
these  techniques  for  presently  elefineel  regions.  Furtlie-rmore, 
analysis  of  another  region  shoulel  be*  implemented  to  asse'ss 
the'  capabilities  of  these  techniques  over  a larger  geograjihical 
event  distribution. 

• When  sufficient  data  are  available,  the  easte'rn  Kazakh  test 
region  should  be'  studied  in  terms  of  these  eiata  e'nhancement 
te'.chniques . 

• SNNR  improvements  for  LQ  shoulel  be  de'termineel. 

• Individual  stations  need  Ui  be  investigateel  in  eletail  in  te'rms 

of  dB  SNNR  improvement  anel  elete'Ction  capability  improve-ment 
due  to  use  of  these  te'chniej  ics. 
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Mean  delay  times  and  associated  standard  deviations  for 
matched  filter  responses  should  be  determined.  These 
are  needed  to  improve  the  detection  criteria. 

Before  the  TCA  processor  is  used  again,  the  phase  and 
true  amplitude  responses  of  the  stations  must  be  determined 
and  corrected  for.  The  question  of  optimum  overlap  and 
gap  length  should  also  be  resolved. 

More  reliance  could  be  put  on  the  detection  levels  if  the 

number  of  test  events  in  the  range  3.  5 < m < 4.  5 were  greater. 

b 

Therefore,  it  is  suggested  that  data  from  a local  network 
or  array  be  used  to  increase  the  number  of  test  events. 

The  PDP-15  interactive  computing  system  should  be 
implemented  to  expedite  the  matched  filter  data  processing. 
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I DETECTION  STATUS 

I In  the  following  tables,  a four-digit  code  gives  the  detection 

status  of  each  station-event  of  the  data  base.  (A  one-digit  zero  code  indicates 
f that  station-event  was  not  included  in  the  data-base.)  The  first  digit  indicates 

the  detection  status  for  the  bandpass  filter,  the  second  digit  indicates  the 
detection  status  for  the  chirp  matched  filter,  the  third  digit  indicates  the 
detection  status  for  the  referenc/^  waveform  matched  filter,  and  the  fourth 
* digit  indicates  the  detection  status  for  the  three-c 'mponent  adaptive  processor. 

A 1 indicates  a detection,  a 2 indicates  a non-detection,  and  a 0 indicates 
I that  particular  data  enhancement  technique  was  not  applied. 
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TABLE  A-1 

CENTRAL  ASIA  DETECTION  STATUS 


EVENT 

NHHBES 

Bfll 

HB4 

885 

886 
887 

895 

896 
900 
902 
905 

908 

909 

910 

911 
913 
Q m 
915 
9 16 

918 

919 

920 
922 
924 
926 
929 
931 
938 
940 
946 
950 

952 

953 

954 

955 
958 
961 
965 
978 

1017 
103  3 
1034 
1044 
1047 
1049 
1051 
1053 

1063 

1064 
1067 
1069 
1084 

1091 

1092 


2 

1111 
1111 
2222 
2222 
21 1 1 
1111 
2212 
2222 
2112 
0 

1111 

1111 

2122 

1111 

0 

1111 

1112 

0 

1111 

2222 
2212 
2112 
2212 
1 r 2 
1111 
2222 
2222 
2222 
1111 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

2110 

2122 

0 

0 

2112 

1112 

1111 


5 

1111 
1111 
2111 
2222 
221  1 
2222 
2222 
2222 
2210 
0 
0 
0 
0 
0 
0 
0 
0 

2122 

1111 

0 

2220 

2120 

2210 

0 

0 

2112 

2212 

2212 

2210 

0 

0 

0 

0 

0 

0 

0 

0 

1111 

0 

2222 

2222 

0 

2222 

2122 

0 

0 

0 

0 

0 

0 

2222 

0 

0 


STATION  NUMBER 
6 8 9 

1111  1111  0 
1111  1111  1111 

2112  2122  2222 
2112  2212  2222 
0 0 0 
1111  0 0 
0 2112  0 
2222  0 0 
222  2 
2122 


0 0 
0 0 
2 0 
0 0 
0 2222 
0 2222 


1111  1111  1111 
2222  2222  0 
2222  2222  2222 
1111  1111  1111 
2222  2112  2212 
2111  2112  2112 
0 2121  2120 
1111  2222  1111 
1111  1111  1111 
0 2222  2222 
0 2222  2122 
O 2222  2222 

0 2122  2222 
1112  2122  0 
1111  1111  0 
2212  2222  0 
2222  212?  0 
2112  2222  0 
2111  1112  0 
nil  nil  1111 
1111  2222  2222 
2112  0 0 
nil  0 1111 

1111  2121  2112 
2122  2222  2222 
2112  2111  2112 
1111  1111  1111 
1111  0 1111 
2111  0 2122 
2111  2212  2212 
2112  2212  2222 
2122  0 2212 


2122 

0 

2222 

2122 

2122 

2221 

2212 

2212 

2122 


0 2222 
0 0 
0 0 
0 2121 
0 2222 
0 0 
0 2222 
0 2222 


2221  2222  2222 
2212  0 2212 
0 0 2212 


10 

1 1 

0 

1111 

0 

1111 

0 

211  2 

0 

2222 

0 

2222 

0 

0 

0 

0 

0 

0 

0 

2222 

2222 

2122 

0 

0 

0 

0 

C 

C 

0 

0 

0 

0 

0 

1111 

2222 

2120 

0 

0 

1111 

1111 

2222 

2122 

2222 

222? 

2222 

0 

2222 

2222 

2222 

1111 

1111 

11  11 

0 

2112 

0 

2212 

0 

2122 

0 

1111 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

nil 

0 

1111 

0 

0 

0 

2112 

0 

2112 

0 

2222 

0 

1111 

0 

0 

0 

0 

0 

2222 

0 

2122 

0 

2112 

0 

2122 

0 

2122 

0 

1111 

0 

2112 

0 

2112 

0 

2222 

TABLE  A-2 

GREECE-TURKEY  DETECTION  STATUS 


EV  KNY 
NOMPER 
RPO 
8^  1 
Rq2 
RR? 
<112 
927 
933 
93U 
939 
9U  1 
94U 
945 
947 

956 

957 
968 

971 

972 
983 

1021 

1028 

1042 

104B 

1052 

1073 

10R1 

10R6 

1087 


2 

1111 

2222 
2212 
211 1 
0 

2121 

0 

222? 

2112 

2222 

2222 

2222 

2222 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2122 

2122 

2222 

2212 


5 

1111 

2221 

2222 

2111 

0 

0 

1110 

n 

2212 

2110 

2220 

2222 

2220 

0 

0 

0 

1111 

0 

1111 

2222 

2222 

2212 

2212 

2222 

2222 

1111 

0 

2222 


STATION  NUMBER 

6 8 9 

1111  1111  1111 
2222  2222  2222 
0 2222  2222 
1111  1111  1111 
1111  0 2212 
1111  2222  0 
1111  2212  0 
0 0 0 
0 0 0 
1111  0 0 


2122  2222  0 
2222  2222  0 
2112  2222  0 

1111  1111  1111 


1111  0 
1111  0 
0 2222 
1111  2222 
1111  0 
0 0 
1111  2122 
1111  0 
0 0 
1111  0 
2112  0 
0 

0 0 
2112  C 


0 0 
0 0 
2 0 
2 0 
0 1111 
0 2222 
2 0 
0 0 
0 2222 
0 2222 
0 2222 
0 2112 
0 2122 
0 2222 


10 

0 

0 

0 

1111 

2212 

2122 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


11 

1111 

2222 

2222 

0 

0 

2112 

0 

2222 

212? 

0 

2222 

2222 

2222 

0 

2222 

0 

0 

0 

0 

0 

0 

2222 

2172 

2222 

0 

2222 
21  1 2 
2122 


APPENDIX  B 


SURFACE-WAVE  MAGNITUDES 


The  following  tables  list  surface -wave  magnitudes  computed 
from  chirp  matched  filter,  reference  waveform  matched  filter,  and  thr^'e 
component  adaptive  processor  data.  These  magnitude  values  are  network 
averages  with  one  or  more  stations  reporting.  A value  of  zero  indicates 
that  no  magnitude  was  computed. 

The  tables  contain  values  only  for  those  station-events 
which  were  detected  b/  the  data  enhancement  techniques  but  not  by  the 
bandpass  filter.  All  matched  fiPor  %aines  were  computed  from  LR-V  data. 
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TABLE  B-1 


FV^’N? 

Rony  wAVF 

S'’PFArF  WAVF 

NTIfIRER 

MAovirnnE 

CUF 

n wi  F 

TCft 

PflS 

4.9 

3.  13 

2.  OP 

3.0  3 

PRR 

3.9 

2.90 

3.10 

0.0 

flB7 

4.7 

2.72 

2.  09 

3.07 

BOB 

4.  .9 

0.0 

2. 87 

0.0 

002 

3.7 

2.  93 

2.  0 3 

0.0 

OOR 

3.8 

3.10 

0.0 

0.9 

01  1 

3.7 

3.  13 

2.  90 

0.0 

0 1U 

4.0 

3.P3 

3.  39 

4.  17 

91B 

4.B 

4.00 

0.0 

5.M 

o?o 

3.7 

2.03 

2.46 

0.0 

922 

4.9 

2.0  2 

2.  82 

0.0 

92« 

4.3 

3.19 

0.9 

0.0 

920 

4.9 

3.23 

C.  9 

0.0 

931 

3.7 

0.0 

3.29 

0.0 

9 3B 

4.0 

3.29 

0.0 

0.0 

940 

9.0 

3.0  9 

O.O 

n.o 

999 

4.2 

3.37 

3.  38 

4.47 

998 

3.9 

2.91 

0.0 

0.0 

901 

4.  3 

3. OB 

3.0  4 

3.  92 

10  17 

4.2 

3.0  3 

3.00 

0.0 

1C13 

4.0 

0.0 

2.  97 

0.0 

1034 

3.7 

3.09 

3.06 

0.0 

1044 

3.4 

0.0 

3.  9 3 

n.c 

1049 

3.0 

2.00 

0.0 

0.0 

1091 

3.6 

2.80 

0.0 

0.0 

10*^3 

9.0 

3.  19 

0.0 

3.69 

1004 

3.R 

4.  19 

4.0  9 

0.0 

1007 

3.9 

2.97 

2.  P3 

0.0 

1 00  9 

3.3 

3.09 

0.0 

9.9 

10P4 

4.9 

0.0 

3.  03 

4.22 

1091 

3.7 

2 . B 2 

3.  28 

0.0 

EVFNT 

pony  WAVE 

SORFAC 

WAVE 

HAO. 

NIIK^FR 

rAONTTHOE 

r^iF 

pw  'IP 

TEA 

897 

S.O 

3.37 

2.  19 

3.11 

912 

9.8 

0.0 

3.  23 

0.  0 

927 

9.8 

3.9p 

0.0 

3.08 

913 

9.7 

o.n 

3.  20 

0.0 

93^ 

9.3 

2.98 

2.  91 

0.0 

991 

9.3 

2.  32 

G.  0 

0.0 

9 9 9 

.9 

2.83 

0.0 

0.0 

997 

3.1 

2.88 

2.  88 

0.0 

1028 

3.8 

3.18 

0.0 

0.0 

1092 

3.7 

0.0 

2.  00 

0.0 

1073 

3.7 

3.06 

0.0 

0.0 

1081 

9.9 

3.  10 

3.81 

0.0 

1088 

9.7 

9.93 

9.  39 

0.0 

1087 

9.0 

2.88 

0.0 

0.0 

